Lipkin's conservation law, Noether's theorem, and the relation to optical helicity 
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A simple conserved quantity for electromagnetic fields in vacuum was discovered by Lipkin in 1964. In recent 
years this "zilch" has been used as a measure of the chirality of light. The conservation of optical zilch is here 
derived from a simple symmetry of the standard electromagnetic action. The symmetry transformation allows 
the identification of circularly polarized plane waves as zilch eigenstates. The same symmetry is present for 
electromagnetism in a homogeneous, dispersive medium, allowing the derivation of the zilch density and flux 
in such a medium. Optical helicity density and flux are also derived for a homogeneous, dispersive medium. 
For monochromatic beams in vacuum, optical zilch is proportional to optical helicity. This monochromatic 
zilch-helicity relation acquires a factor of the square of the phase index in a dispersive medium. 



PACS numbers: 03.50.De, 42.5().Wk 



I. INTRODUCTION 



Field theories such as electromagnetism possess an infinite 
number of conserved quantities in the absence of sources iHJ- 
@]. One can view electromagnetic conservation laws as con- 
straints on the field dynamics and, in this sense, every con- 
served quantity is physically significant: the list of conserved 
quantities is an alternative statement of information contained 
in Maxwell's equations. On the other hand, only the informa- 
tion contained in a limited number of conserved quantities has 
a straightforward relation to basic physical concepts. In 1964 
Lipkin [1] discovered an electromagnetic conserved quantity 
that is as simple as those related to the basic concepts of en- 
ergy, momentum and angular momentum. The density and 
flux of Lipkin's "zilch", as he termed it, are 



Px =f2 (-E -d t B + B. 3 t E) , 
S x =^-E x d t E + — !— .B x d t B. 
These obey the conservation law 



0, 



tPx 



s x = o, 



(1) 

(2) 
(3) 



showing that the integral of the density p x over all space is 
conserved. In recent years p x has been used as a measure 
of the chirality density of light, giving a simple method of 
exploring light configurations that couple to material chiral- 
ity I7H13II. One example of a light beam with a non-zero zilch 
density p x is a circularly polarized plane wave, and this has 
motivated researchers to compare and contrast zilch with the 
spin angula r momentum of light and the helicity of the pho- 
ton j 14141911 . It is important to note that there are more in- 
teresting possibilities than circular polarization for non-zero 
values of p x , once local variations are considered. In par- 
ticular, it has been shown that the evanescent fields of metal 
patterns can be engineered to produce an enhanced local zilch 
density Q [Til UM and light beams with orbital angular mo- 
mentum can have a locally enhanced p x on the axis of the 
beam 1 1 311 - 



The conservation of zilch must be associated with a sym- 
metry of the electromagnetic action 

S[A,4>] = ^ I ' d i x{E.E-c 2 B.B), (4) 

which is a functional of the vector potential A(r, t) and the 
scalar potential <p(r,t) with 



d t A, B = V X A. 



(5) 
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Once a symmetry of an action is found, Noether's theorem 
provides a constructive proof of the associated conservation 
law (see 12011 . for example). If however a conservation law 
is discerned using the field equations, as was the case with 
Lipkin's discovery QJJ] , then the identification of the associ- 
ated symmetry of the action may not be straightforward. The 
relation of zilch to a symmetry transformation has been con- 
sidered several times before IEl1 - l28l . but apart from [5TI the 
transformations considered were not of the scalar and vec- 
tor potential in the standard electromagnetic action (@). In- 
stead actions with additional dynamical variables were con- 
sidered I22I - I28I1 . for example functionals of the vector poten- 
tial, scalar potential and the electric field, functionals of the 
transverse vector potential and the electric field, or function- 
als of the vector potential and an "electric" vector potential. 
In some cases the particular non-standard action is solely mo- 
tivated by the desire to derive the conservation law from an 
action as easily as possible j22l - l26Tl . In other cases the form 
of the non-standard action is motivated by an interest in deriv- 
ing a conservation law associated with electric-magnetic du- 
ality I27I.I28I1 . the conserved quantity in this case being optical 
helicity. Incorporating electric-magnetic duality into the stan- 
dard action presents some difficulties II27[ 12811 . as can be seen 
in [22] where the relevant symmetry transformation contains 
an inverse Laplacian operator (and the action was also rewrit- 
ten as a functional of the transverse vector potential and the 
electric field). Calkin's brief note [21] does not directly de- 
rive zilch conservation from a transformation of the dynami- 
cal variables A and <f> of the standard action (0), as discussed 
in Sec mi 

Our goal here is to derive zilch conservation from a sym- 
metry of the standard action (|4) using Noether's theorem. As 
the action © is a functional of the independent degrees of 
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freedom of the electromagnetic field up to gauge invariance, 
the symmetry underlying zilch conservation takes its simplest 
form in terms of these dynamical variables (Sec.|Il). The sym- 
metry is used in Sec. [HI] to construct zilch eigenstates, in a 
manner analogous to energy, momentum and angular momen- 
tum eigenstates. Section [IV] outlines generalizations of the 
symmetry transformation that lead to an infinite class of con- 
served quantities |Q]-§]. In Sec.lVlwe consider the electromag- 
netic action in a homogeneous, dispersive medium; the zilch 
symmetry is present for this action and we use it to find the 
zilch density and flux in such a medium. The relation of zilch 
and helicity to the notion of optical chirality is still debated 0- 
H^l . In Sec.[VT]we generalize the expressions for optical helic- 
ity density and flux to a homogeneous, dispersive medium and 
show that zilch, helicity and spin angular momentum all have 
a different dependence on the material dielectric functions. 



II. OPTICAL ZILCH IN VACUUM 



It is easy to show that the vacuum electromagnetic action 
© is invariant (up to surface terms) under the following active 
transformation of the dynamical variables A(r,t) and </>(r, t): 



6A = 77V x d t A, 



= 0, 



(6) 



where 77 is the infinitesimal parameter of the transformation. 
Two unusual features of the transformation © compared to 
the Poincare-group symmetries of © are the absence of any 
change in <fi and the occurrence of two derivatives in the in- 
finitesimal transformation of A. If either of the derivatives in 
5 A are removed the result is not a symmetry transformation, 
but any even number of additional space and time derivatives 
can be inserted to give new symmetry transformations (see 
Sec. |IV| ). Calkin 112111 considered transformations of the form 
5 A = r\ V X Z for any Z satisfying V • Z = 0. Note that such 
a transformation does not include © because the three com- 
ponents of A are independent in the action ©. Calkin's ap- 
proach allows a quick extraction of optical zilch from the ac- 
tion ©, but here we seek the conservation law from a symme- 
try transformation of the dynamical variables of the action, as 
is done for energy-momentum and angular momentum ll20ll . 

Noether's theorem [20] shows that if we let the parameter r\ 
in the symmetry transformation © depend on space and time, 
i.e. r\ = ?7(r, t), then the change in the action © under © can 
be written in the form 



5S = J d 4 x {p x d t T] + S x • Vrj) , 



(7) 



where p x and S x obey the conservation law ©. We now 
show that p x and S x in © are optical zilch density and flux, 
respectively. After various integrations by parts (surface terms 



are dropped) the form © is achieved with 
P X =eo (-E ■ d t B - ^d t B ■ d t A^j - ^-A • V 2 B, 



(8) 



S x =e [ -4>dfB + -dtA X Of A 



1 

Mo 



d t B X B - (V • d t A)B - t; a x v d * A 



(9) 



We use tensor notation and the summation convention, so that 
all index summations are contractions of an upper and a lower 
index. One can verify that © and © satisfy the conservation 
law © because of the dynamical equations of the action (01), 
namely Maxwell's equations 



(10) 



(The other two Maxwell equations are identities because of 
©•) 

The conserved quantities that emerge directly from 
Noether's theorem applied to the action © are not gauge in- 
variant 12011 , Gauge-invariant densities and fluxes are found 
by the addition of terms that identically satisfy the conser- 
vation law. As regards densities such as p x , the divergence 
of any field quantity X(A, <fi) can be added without chang- 
ing the value of the conserved quantity J d 3 r p x ; if the time 
derivative of X(A, (f) is then subtracted from S x , the conser- 
vation law (f3]) is preserved. One verifies that the addition to 
© of the divergence 



■^.tpdtB* + —AX 'IS* 
2 p 



(11) 



where the square brackets denote anti-symmetrization (i.e. 
TM = (TV - T ji )/2), produces the zilch density ©. To 
preserve the conservation law the time derivative of the quan- 
tity in brackets in ( fTTT i must be subtracted from ©; this does 
not produce a gauge-invariant flux but any term with zero di- 
vergence can now be added to S x while maintaining ©. The 
simplest course is to show that, after the subtraction just de- 
scribed, the divergence of © is equal to the divergence of © 
(without use of Maxwell's equations (fTUl)). This demonstrates 
that © and © are gauge-invariant versions of © and ©, 
respectively. 

The symmetry transformation © is simpler than transfor- 
mations that have previously been associated with the con- 
servation of optical zilch Il2ll - l28ll . This is because © is a 
symmetry of the standard electromagnetic action © that is a 
functional of the independent degrees of freedom of electro- 
magnetism up to gauge transformations. The standard quanti- 
zation of electromagnetism is based on © and it is most natu- 
ral to relate conservation of optical zilch to a symmetry of this 
action, as is done for conservation of energy-momentum and 
angular momentum [20]. 
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III. OPTICAL ZILCH EIGENSTATES 

The infinitesimal versions of continuous symmetry trans- 
formations reveal the generators of the symmetry group asso- 
ciated with a conserved quantity. Eigenstates of these gener- 
ators have eigenvalues that completely determine the value of 
the conserved quantity carried by the eigenstate. This is most 
familiar in quantum mechanics but the same procedure can be 
applied to classical field theories. The generators {idt, iV} of 
space-time translation symmetry give plane waves e lk " r - lut 
as energy-momentum eigenstates. The components of the 
vector generator L = ir X V of rotations do not commute 
but eigenstates of L 2 can be found and these give angular- 
momentum eigenstates l29ll . The eigenstates are taken to be 
complex, with the electromagnetic fields given by the real part 
of the eigenstates. From (|6]l the generator of the symmetry as- 
sociated with optical zilch is idt V X and it acts on the vector 
potential. We find the vector-potential zilch eigenstates by 
considering the Fourier-transformed field A(r, ui) in the fre- 
quency domain; note that this field is complex for real time- 
domain fields A(r, t), so we do not need to consider complex 
A(r, t) to find the eigenstates. The zilch eigenvalue equation 
for A(r, to) is then 

V X A(r,u) = ^A(r,u), (12) 

UJ 

where \ is the eigenvalue. We consider a gauge in which <j> = 
0. Then 

E(r,u) =iu)A(r,uj) (13) 

and ( fT2l gives 

B(r,w) = -i-^E(r,u>). (14) 

UJ 1 

Equation ( flZb states directly that B(r,u) = (x/w)A(r,uj); 
taking the curl of this and using (fT2l gives 

2 

V X B(r,uj) = X A(r,u) = ^l(r,w) 

u> U! z 

= -i*-E(r,u), (15) 

where the last step uses ( fT3l . Comparing (fTBI l with the second 
Maxwell equation ( fTOb in the frequency domain, we see that 
the eigenvalue \ must be 

X = ± — ■ (16) 

c 

Insertion of ([Tol l in (fl4l gives 

B(r,u>) = T~E(r,u>), (17) 

c 

which is precisely the relation between the complex ampli- 
tudes of the electric and magnetic fields in a circularly polar- 
ized plane wave [29]. We did not assume that the zilch eigen- 
states are monochromatic — no restrictions were placed on the 



frequency-domain spectrum A(r, w). We see from (fTST l how- 
ever that x depends on frequency, so for it to be an eigenvalue 
of the total field in the time domain we can only have one fre- 
quency component. Moreover (fTTI l shows that this monochro- 
matic wave must be circularly polarized. Circularly polar- 
ized plane waves are thus zilch eigenstates with eigenvalues 
JTol l. Note that there is degeneracy in the eigenstates: only the 
frequency, not the direction, of the circularly polarized plane 
wave determines the eigenvalue ( fTol l. Any superposition of 
plane waves with the same frequency and circular polariza- 
tion is also a zilch eigenstate with eigenvalue ( fTot . 

If we replace the zilch generator with its quantum version, 
namely ihdt'Vx, the eigenvalues ( fTBT l become 

X = ±^- (18) 

c 

The eigenvalues ( fT8l per photon energy hoj are then ±ui /c, 
which is well known to be the time-averaged zilch per unit 
energy of a classical circularly polarized plane wave 1171 41311 . 

It is not surprising that the zilch eigenstates should dis- 
tribute the zilch density uniformly throughout the wave, which 
is what the circularly polarized plane wave does. For appli- 
cations however, beams with nonuniform, locally enhanced 
zilch density are more interesting ||8l llll — Tl3l1 . 

In Sec. [V] the vacuum zilch eigenvalues ( fT6] > and the vac- 
uum time-averaged results for the zilch of plane waves will be 
generalized to the case of a dispersive medium. 



IV. GENERALIZATIONS OF THE SYMMETRY 
TRANSFORMATION 

As noted in Sec. [TTJ the transformation ((6) is still a symme- 
try of the action when any even number of space and time 
derivatives are inserted into the transformation of A. The total 
number of extra derivatives inserted must be even, so an odd 
number of space derivatives together with an odd number of 
time derivatives is allowed. Depending on the number of extra 
space derivatives and whether the associated indices have any 
contractions, we can construct an infinite number of conserved 
quantities with any number of tensor indices. The resulting 
infinite class of conservation laws were found by Morgan (U. 
As an example we note that the generalized symmetries that 
give a infinite class of scalar conserved quantities are 

5A = V V x V 2n 8 2m A, 5<t> - 0, (19) 

where n and m are any non-negative integers. The conserva- 
tion law associated with ( fT9T > is calculated as in the n = m = 
case. After some tedious manipulations to bring the density 
and flux into their simplest gauge-invariant forms, we find the 
conservation law 

9tP(n,m) + V • S( n ,m) = 0, (20) 
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with 



(21) 



Si n m ) =^-E x W 2n d 2m+1 E + ^-B x V 2n d 2 t m+1 B. 
2 2,uo 



(22) 



It is straightforward to verify that (fJT) and (l22l satisfy (l20l > 
with use of Maxwell's equations dTOb . The Laplacian oper- 
ators V 2 in (fJTJ and (l22l can be replaced by the operator 
(l/c 2 )<9 2 because of the wave equation. 



V. OPTICAL ZILCH IN A DISPERSIVE MEDIUM 

Electromagnetic waves in materials can carry conserved 
quantities in limited frequency ranges where absorption is 
negligible. The effects of dispersion however can be impor- 
tant even when absorption can be ignored. In this regime of 
negligible absorption the best known example of a conserved 
quantity in a dispersive medium is the energy of a monochro- 
matic wave, whose time-averaged energy density is given by 
Brillouin's formula (l38l below (see 112911 . for example). The 
electromagnetic energy density has a more complicated ex- 
pression for waves that are not monochromatic [30]. The 
electromagnetic momentum density [30] and angular momen- 
tum density [31] that measure the conserved values of these 
quantities inside the medium also have dispersive contribu- 
tions, which take a simple form in the monochromatic case. 
The conserved electromagnetic energy-momentum and angu- 
lar momentum in a dispersive medium can be calculated using 
Noether's theorem BOiblll . and we now do the same for elec- 
tromagnetic zilch. 

The dielectric functions in a limited frequency range where 
absorption is negligible can be fitted to an even series in fre- 
quency (see 113211 . for example): 



e(w) = ^2 e 2n o 



,2n 



71=0 



E 

ra=0 



.-a 



(23) 



where e(u) is the relative permittivity and the relative perme- 
ability is = k(cj) -1 . The series in (l23l are a fit to the 
dispersion data of the material in the frequency range of inter- 
est and as such will be finite; we do not place an upper limit 
on the summations as the results will be valid for infinite se- 
ries. We assume the medium is homogeneous as this will be a 
requirement for conservation of optical zilch. The D and H 
fields in the frequency and time domains are 

D(r,w) = e e(w)^(r,o;), D(r,t) = e e(id t )E(r,t), 
H{r,u>) = K K(uj)B(r,uj), H(r,t) = K K(id t )B(r,t), 

where kq = fj,^ . The electromagnetic action in the homoge- 
neous dispersive medium ( 1231 is j30l[3lll 



S[A,< 



d 4 x 



K f 1 



-E ■ [e(id t )E] - B ■ [K{id t )B] 



Variation of and A in d24l gives the macroscopic Maxwell 
equations 



e V ■ [e(id t )E] = 0, 
k V x [K(idt)B] = s s{id t )d t E. 



(25) 
(26) 

The action (1241 is still symmetric under the transformation 
© (this transformation is not a symmetry if the medium is 
inhomogeneous, i.e. if s(idt) e(r,idt) and/or n(idt) —> 
K(r,id t )). The conservation law associated with the sym- 
metry (O is extracted from the action (l24l using Noether's 
theorem, as in the vacuum case of Sec. [XT] There is now the 
extra complication of the series in (l24l . The method of deal- 
ing with these series when applying Noether's theorem is de- 
scribed in IBOIl and 03111 . to which we refer the reader for the 
details. The zilch density and flux in the resulting conserva- 
tion law can be made gauge invariant as in Sec. [II] and the 
result is 

Px=-\d- d t B + l -B ■ d t D 

oo 2n 

+f E E (-ir +m £ 2 „sr£ . a* n - m v x e 

n— 1 m—1 

oo In 

+ Y E E (-l)" +m K 2 n9 t m " 1 -B • d 2n - m+l V X B, 
n—1 m—1 

(27) 
(28) 



S y =-D x d t E + -H x d t B 

x 2 2 



It is straightforward to show that (l27l and (128} obey the con- 
servation law (0 when the macroscopic Maxwell equations 
(|25l and d26l are used. 

As in the cases of energy, momentum and angular momen- 
tum POL 13111 . the expression for the zilch density (l27l in a dis- 
persive medium simplifies considerably for a time-averaged 
monochromatic beam. We consider the monochromatic E 
field 



E{r,t) = ±(E (r)e- i " t + C . C ) 



(29) 



and a B field of the same form. Inserting these into (1271 and 
(|28l i and taking a time average we obtain the monochromatic 
time-averaged zilch density p x and flux S x : 

E s(ui)uJ 

Px = s MJEo-Bo) 



+ 



2 

Equj de(uj) 

4 doj 
KqUJ dK(uj) 

4 du 



Re{E ■ V X E*) 
Re(B .VxB *), (30) 



Im(B X B*) 



(31) 



(24) 



The dispersive contributions to the time-averaged monochro- 
matic zilch density d30l contain first derivatives of the dielec- 
tric functions at the frequency of the beam; similar dispersive 
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contributions occur in the momentum density rf30ll . angular 
momentum density 03111 . and energy density (Eq. ( 138) below). 
The curls in d30b can be rewritten using the Maxwell equa- 
tions V X E = iuB , k(w)V X B = —iuje(u})E a /c 2 for 
the monochromatic complex amplitudes. Utilizing also the 
definitions 



d[ujn p (ui)} 
dui 



n p (u) = y/E^jjfMuj), (32) 



of the group index n g (uj) and the phase index n p (uj) at the 
frequency of the wave, we find a simple formula for the time- 
averaged monochromatic zilch density (130) : 



The zilch per unit energy of a circularly polarized plane wave 
is thus proportional to the phase index of the medium. This 
is in sharp contrast to the conserved spin angular momentum 
per unit energy of a circularly polarized monochromatic beam, 
which is independent of material properties even in a disper- 
sive medium j3lll . On the other hand, the conserved linear 
momentum per unit energy of a plane wave is also propor- 
tional to the refractive index in a dispersive medium Dill . 

The zilch eigenstates of Sec. [HI] are easily generalized to a 
dispersive medium; circularly polarized plane waves are again 
eigenstates but the eigenvalues now acquire a factor of the re- 
fractive index: 



(33) 



°x — 



s u 



e(u)lm(E X E*) 



Tm (Bq X B*) 



(34) 



where we restate the time-averaged monochromatic zilch flux 

(ED. 

We apply the monochromatic results to a circularly polar- 
ized plane wave propagating in the x direction. The complex 
amplitudes Eq and Bq of the electric and magnetic fields in 
such a beam are 1 29] 



E (r) = £(e y ± ie z )e 



k = -njoj), 
c 



ikx 



B (r) = T-n p (u)E (r), 

c 

(35) 



n p (u) = x/e(uj)p(u), 



(36) 



where e y (e z ) are unit vectors in the y (z) directions, k is the 
wave number, and n p (ui) is the phase index as in d32) . The 
time-averaged zilch density and flux of this wave are, from 
(03) and (03), 

Px = ± l ^e e(Lu)n g (uj)£ 2 , S x = ±Lue e{uj)£ 2 e x . (37) 

Dividing S x by p x in d37) we find that the optical zilch flows 
through the medium at the group velocity c/n g (u)), just like 
the optical energy. To find the time-averaged zilch per unit 
energy in the circularly polarized plane wave, we require the 
time-averaged energy density p of a monochromatic beam in 
a dispersive medium, which is given by Brillouin's expres- 
sion 129, 3JJl 



£o d[uje(uj)} 



For the plane wave 
i is 



E* 



Ho d[u>fj,(u)} 



H H*. (38) 



4 duj 

, the time-averaged energy density 



e n p (uj)n g (LL>) 
p(uj) 



(39) 



Dividing p x in ((37) by d39b we find the zilch per unit energy: 

±k. (40) 



Px 
P 



±^-n p (uj) 



X = ±—n p (uj). 



(41) 



VI. OPTICAL HELICITY IN A DISPERSIVE MEDIUM 

A topic of recent discussion is whether the notion of op- 
tical chirality is best defined through zilch or through helic- 
ity lfl4l - [l9ll . While zilch and helicity are different quantities 
with different dimensions, they are proportional to each other 
for monochromatic beams in vacuum (and only for monochro- 
matic beams). It is interesting to extend the comparison of 
zilch and helicity to a homogeneous, dispersive medium. This 
requires the expressions for optical helicity density and flux in 
the dispersive regime considered in the previous section. 

Optical helicity has its natural expression in terms of mag- 
netic and electric vector potentials 11141 12211 . As long as the 
medium is homogeneous we have V • E(r, t) = and so we 
can introduce an "electric" vector potential C(r, t) as follows: 



E 



-cV x C. 



(42) 



Maxwell's equations in the regime considered in the previous 
section then give the following relations between the electric 
(magnetic) field and the magnetic (electric) vector potential: 

E = -d t A, n(id t )B = --e{idt)d t C, (43) 

c 

We first note that it is easy to generalize optical helicity den- 
sity and flux to a non-dispersive medium where the dielectric 
functions are constants (rather than operators) in the time do- 
main. The expressions are obtained from the vacuum helicity 
density and flux iffll^ by inserting appropriate factors of the 
non-dispersive e and k(= l/p) (here we use SI units through- 
out): 



Ph 



£0 



{ckA -B- eC-E) 



Sh =^Y K (Ex A + cB xC). 



(44) 
(45) 



These expression were used in 113 311 for a monochromatic 
beam. It is easy to show that (04) and (05) obey the con- 
servation law 



dtp h + V • S h = 



(46) 
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because of the non-dispersive Maxwell equations (and rela- 
tion between E, B and the vector potentials). On the other 
hand, it is also easy to see that the conservation law (l46l l is 
not satisfied by (l44b and d43T > if the medium is dispersive be- 
cause e and k are then operators and the dispersive Maxwell 
equations do not give d46b no matter which vector fields the 
operators e{idt) and n(idt) are taken to act on in ( |44T > and 
d45l l. The situation here is no different than in the familiar 
case of electromagnetic energy 1I29I l30ll : the energy density 
in a non-dispersive medium does not generalize in a simple 
way to a dispersive medium. There are non-trivial dispersive 
contributions to electromagnetic energy, momentum and an- 
gular momentum in media II29I431I1 : we have seen that the 
same is true for zilch and now we must supply the disper- 
sive contributions to helicity. The systematic way of finding 
the correct dispersive contributions is through Noether's the- 
orem, as in Sec. [V] In the case of helicity however, the use 
of the standard electromagnetic action (@) presents some dif- 
ficulties 1221 l27ll28ll . Experience gained through solving the 
cases of energy-momentum, angular momentum and zilch al- 
lows us to simply write down the correct dispersive contribu- 
tions to the helicity density. Whereas finding the generaliza- 
tion of a electromagnetic conserved quantity to a dispersive 
medium can be challenging, it is a straightforward matter to 
check whether any given generalization is correct, i.e. whether 
it satisfies the conservation law d46| i. Consider the expressions 



Ph =f|cA • [K(id t )B] - [e(id t )C] ■ E 

00 In 

+ ^2 Y (-i)™ +m £2nar _1 c • d 2 t n - m+i E 

n=l m=l 

00 2n ^ 

+c Y Y (-l)" +m K 2 „a t m " 1 A • d 2 t n ~ m+1 B \ , (47) 

n— 1 m— 1 ) 

S h = °^- {E X [n{id t )A] + c [K(id t )B] X C} . (48) 



These satisfy d46b with use of the macroscopic Maxwell equa- 
tions and the relations (l42l and (@3j- The expressions ( |47t 
and d48l i also give vacuum helicity density and flux as a spe- 
cial case; they are therefore the correct generalization of con- 
served optical helicity to a homogeneous dispersive medium. 



For a monochromatic wave (129) , the time-averaged helic- 
ity density and flux are easily found from (l47l i and (|48V The 
calculation is almost identical to that for zilch and we aid the 
comparison by eliminating the complex amplitudes of the vec- 
tor potentials using 



lo — E a , Cq 

U! 



-1 — rr^Oi 

(jJE(LU) 



(49) 



averaged helicity density ph and flux Sh as 

ce n g (uj) 



Ph 
S h 



2uj p{u>)n p {ui) 



lm(E Q -B* 



0) > 



(50) 



e(uj)lm(E X E*) 



p(uj) 



lm(B X B*] 



(51) 



Comparing these with the corresponding results (133) and (l34l 
for zilch, we find the following proportionality between helic- 
ity and zilch for any time-averaged monochromatic wave in a 
dispersive medium: 



Ph 



cj 2 n2 (uj) 



Px> 



aj 2 nl(uj) x ' 



(52) 



This proportionality is familiar in the vacuum case n p (lu) = 1; 
we have shown that it contains the square of the phase index 
in a dispersive medium. 

The zilch results (|37| | for the circularly polarized plane 
wave (135) . together with the proportionality (l52t . immediately 
give us the time-averaged helicity density and flux of a circu- 
larly polarized plane wave in a dispersive medium: 



s a n g {uj) c2 



S-y i . . 

Ojp(UJ) 



(53) 



The proportionality (1521 implies that, like zilch, helicity 
moves through the medium at the group velocity. The helicity 
per unit energy is given by p x in (l53l divided by 1 



Ph 
P 



±- 



1 



ujn p (uj) 



(54) 



which follow from (143 1 . This gives the monochromatic time- 



The helicity per unit energy of a circularly polarized plane 
wave is thus inversely proportional to the phase index n p (w), 
whereas the zilch per unit energy (l40l is proportional to 
n p {uj). Although helicity has the dimensions of angular mo- 
mentum, the (spin) angular momentum per unit energy of a 
circularly polarized monochromatic beam is ±1/cj even in a 
dispersive medium i3~Ill and so behaves quite differently to 
helicity per unit energy (|54l . 



VII. CONCLUSIONS 

The symmetry underlying conservation of Lipkin's zilch is 
a simple transformation of the vector potential (Eq. (|6j). Iden- 
tification of the symmetry transformation has allowed us to 
prove that circularly polarized plane waves are optical zilch 
eigenstates. These plane-wave eigenstates, however, are not 
the most interestinepossibility for light beams that carry op- 
tical zilch J^, [TTl - [l3ll . A straightforward generalization of 
the symmetry transformation gives the symmetries underly- 
ing an infinite class of electromagnetic conserved quantities 
identified by Morgan [2]. Optical zilch and optical helicity 
are both conserved in a homogeneous, dispersive medium, in 
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frequency ranges where absorption is negligible. For time- 
averaged monochromatic waves, zilch and helicity differ by a 
factor that contains the square of the phase index of the ma- 
terial. The time-averaged zilch per unit energy of a circularly 
polarized plane wave in a dispersive medium is proportional 
to the phase index, whereas the helicity per unit energy is in- 
versely proportional to the phase index. In contrast to zilch 
and helicity, the conserved (spin) angular momentum per unit 
energy of a circularly polarized monochromatic beam in a 
dispersive medium is independent of the dielectric functions. 



These results demonstrate further contrasts between zilch, he- 
licity, and spin angular momentum. 
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